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INTRODUCTION 



k energy problem, which is causing the accelerating expansion of the Universe in the large 
The most accepted observational evidence for existence of such component comes 



As we know Dark (non-luminous and non-absorbing) Matter (DM) is an old idea even before 
the darl 
scale 

from the astrophysical measurements of several galactic rotation curves. From the view of classical 
mechanics, we expect that the rotational velocity of any astrophysical object moving in a stable 
(quasi stable) Newtonian circular orbit with radius r must be in the form 



%(r) oc v / M(r)/r, 

where M(r) is identified as the mass (effective mass) profile thoroughly inside the orbit. For many 
spiral and elliptical galaxies this velocity v v remains approximately constant for large galactic 
radii, for instance in the Milky way galaxy, v ~ 240km/ s. This estimation is valid only near 
the position of our solar system. There is a lower bound on the DM mass density, Odm ~ 0.1 
from phenomenological particle physics. There are different kinds of the dark matter To solve 
the DM problem several proposals were introduced. The problem can be interpreted as an effect 
of the extra dimensions in a cosmological special relativity (CSR) model, proposed by Carmelli 



From particle physics view as LSP in super symmetric theories or LKP in higher dimensional 
theories in which the SM (standard model) predicts some extra dimensions. Stability condition 
of any candidate for dark matter is very important problem which must be checked. For example 
for stabilization checking in SUSY (super symmetry) we must check the validity of R- parity and 
in supergravities alternatives we must following the KK parity. In brief, 11 Any candidate for dark 
matter need not be stable if its abundance at the time of its decay is su fficiently small" . There 
are several classical candidates for dark matter, as perfect fluid models [5[ and as the geometrical 
modifications of the Einstein-Hilbert action, for example R 2 modification of the usual Einstein 
gravity [g] or in the anisotropic and dfiffeomorphism invariance model of Horava-Lifshitz as an 
integration constant Q|. 

In this Letter we focus on the mechanism of the f(T) gravity and show that in the context of 
this new proposed non-Riemannian extension of the general relativity (GR), it is possible to explain 
the rotation curves of the galaxies without introducing dark matter. Our plan in this letter is: In 
section II we propose the basis of the f(T) gravity. In section III, we investigate the spherically 
symmetric solutions of the model. In section IV we solve the equations and show that the rotation 
curve of the galaxies in this toy model of the spherically-symmetric-static model can be recovered 



3 



by the effects of the torsion alone. In section V we obtain the halo density profile and compare it 
with two well known astrophysical models. We conclude in the final section. 

II. FORMALISM OF f(T) GRAVITY 

A gauge theory of gravity is based on the equivalence principle. For example, SL(2 1 C) gauge 
theory on the gravitational field can be used for quantization of this fundamental force [8(]. We are 
working with a curved manifold for the construction of a gauge theory for gravitational field. It is 
not necessary to use only the Riemannian manifolds. The general form of a gauge theory for gravity, 
with metric, non-metricity and torsion can be constructed easily [9]. If we relax the non-metricity, 
our theory is defined on Weitzenbock spacetime, with torsion and with zero local Riemann tensor 
R^fa = 0. In this theory, which is called teleparallel gravity, we use a non-Riemannian spacetime 
manifold. The dynamics of the metric determined using the torsion T. The basic quantities in 
teleparallel or the natural extension of it, namely f(T) gravity is the vierbein (tetrad) basis 



1CH12I]. This basis is an orthonormal, coordinate free basis, defined by the following equation 



9pv — e p e i r lij 



This tetrad basis must be orthonormal and rjij is the Minkowski flat tensor. It means that e^ e^ = Sj. 
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One suitable form of the action for f(T) gravity in Weitzenbock spacetime is given by 

S = j d*xe(^[T + f(T)} + L m ), (1) 
where / is an arbitrary function, e = det(e!J. Here T is defined by 

rp o pVrpp 

1 ^ p 1 pv 1 

with 

T P pv = e f(^ e ^ — ®v e p) > 

Sp ,» = 1 _ {K ,v p + 6 » T e»_ 6 v T ep h 
where the asymmetric tensor (which is also called the contorsion tensor) K^ 1 ' p reads 
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The equation of motion derived from the action, by varying the action with respect to the e*, is 
given by 

e-'d^eS, H(l + fir) ~ ^S/" fir + S t ^d,(T)f TT - \e\(\ + /(T)) = AnefT/ (2) 

where T^ u = e^T a „ is the energy- momentum tensor for matter sector of the Lagrangian L m , it is 
defined using 

m 1 5L m 



e 5e Ql/ 

The covariant derivatives compatible with the metricity g-^ = 0. It is a straightforward calculation 
to show that ([2]) is reduced to the Einstein gravit y w hen f(T) = 0. This is the equivalence between 



the teleparallel theory and the Einstein gravity 



1 31 ] . Note that teleparallel gravity is not unique, 



since it can either be described by any Lagrangian which remains invariance under the local or 
global Lorentz SO(3, 1) group [uj]. 

We mention here that a general Poincare gauge invariance model for gravity (is so called 



Einstein-Cartan-Sciama-Kibble (ECSK)) previously reported in the literatures 15(. Specially, in 
the framework of the Poincar gauge invariant form of the ECSK, theory the notions of "dark mat- 
ter" and "dark energy" play the role similar to that of "ether" in physics before the creation of 
special relativity theory by Einstein [16]. In this Letter we focus only on f(T) models, without 
curvature and with non zero torsion. We would remark that an attempt to explain the flat 
rot ation curve of galaxies has been made earlier in the framework of ECSK theory 



17] 



III. SPHERICALLY SYMMETRIC GEOMETRY 



As is well-known that a typical spiral galaxy contains two forms of matter: luminous matter 
in the form of stars and stellar clusters which are found in the galactic disk while another is dark 
matter which is generally found in the galactic halo and encapsulates the galaxy disk. In the early 
Universe, the DM played crucial role in the formation of galaxies when the dense concentration of 
DM lied in the galactic centers which helped in the accumulation of more dust and gas to form 
proto-galaxies. In the later stages of galactic evolution, the DM slowly drifted towards the outer 
regions of the galaxies forming a huge (but less concentrated/dense) DM halos. Although the 
precise form of distribution of dark matter in the halos is not known, we assume that the spatial 
geometry of galactic halo is spherically symmetric. Moreover, the dark matter halo is isotropic: the 
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spherical DM halo expands (hypothetically) only radially while having no tangential or orthogonal 
motions relative to the radial one. From the point of view of Grand Unified Theories (GUT), the 

n 

most likely candidate of DM is neutralino which is a weakly interacting massive particle |18|] with 
additional minor contribution form primordial black holes S^pbh = 1CU 8 which were formed in the 
early Universe and are also candidate for other violent cosmic events like Gamma Ray Bursts [it] ]. 
Note that we are not interested in any particular form of DM and deal only with its characteristic 
role in the rotation of galactic disks. Our theoretical model suggests that the flat galactic rotation 
curves can be explained in terms of torsion of space without invoking dark matter. In other words, 
the huge DM halo is nothing but mysterious and elusive torsion of space. Now we construct a 
model for galaxy, based on the above assumptions. The metric of a static spherically symmetric 
(SSS) spacetimes can be described, without loss of generality, as 



ds 2 = e a ^ dt 2 - dr 2 - r 2 (d6 2 + sin 2 0# 2 ) . (3) 

This form of metric is updated by the Schwarzschild gauge, and is useful to construction of a toy 
model for galaxy. In order to re-write the metric ([3|) into the invariant form under the Lorentz 
transformations, we use the tetrad matrix j^cl ] 

{e^} =diag{e a(r)/2 ,e b{r)/2 ,r,rsme} . (4) 

Although the f(T) gravity is not local Lorentz invariance [2l|], but we can impose local invariance 
symmetry on the metric components. Here we note some remarks about the relation between local 
Lorentz invariance of the f(T) as a scalar gravitational theory and choice of the tetrads, specially 
in the case of spherically symmetric metric, given by (j4|). As we know, teleparallel gravity with 
f(T) = T is local lorentz invariance, for any set of the constants C12, since it is equivalence to 
the Einstein theory, it means the total Lagrangian of the linear T theory is equivalence to the 
Einstein pi 
of motion 

free of the dynamics. Torsion based f(T) theory, constructed from the tetrad basis must be 
local invariance under a proper Lorentz transformations. But in this form and with usual tetrad 



us a surface boundary term which can be canceled in the derivation of the equations 
ljj. Boundary terms like n u Tff u (T is torsion) have thermodynamical meaning but are 



211 ] . Recently (without any direct proof), the 



basis, it has been shown that, this invariance breaks 
authors of 22], shown that with another choice of the tetrads (they called "good tetrad"), which 



leads to the non diagonal metric components, the restriction on the form of the f(T) as a linear T 
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theory, came from the equation 



frr = 0, 



relaxes. This relaxation can be interpreted as a rotation of the tetrads basis. In this new tetrads 
basis, there are three free Euler angels (all local functions of the coordinates r, t), and the expression 
of the scalar torsion has some extra terms. In this case the system of equations even for vacuum case 
is very complicated. Further, they found, it is possible to recover the Schwarzschild-de Sitter as an 
exact black hole solution in this theory, and test the parameters of the model, using the usual tests 

n, 

as PPN and etc. There are some points that must be clarified by those authors [22| to support their 
conclusions, for example does this local invariance help to the power counting renormalizability of 
the f(T) model as an alternative theory? There are many kinds of such rotational transformations. 
How we preferred one kind of the rotations from others ? We will be back to these problems in a 
forthcoming paper on this topic [23]. Any way in this letter we want to explore the effect of the 
torsion for construction of a geometrical model for DM. For this purpose, the usual tetrads are 
enough. 

Using one can obtain e = det [e*„] = e ( a + b )/ 2 r 2 s i n (#), and the torsion scalar in terms of r 
is given by 



T(r) 



2e" 



a + 



1 



(5) 



where the prime (') denote the derivative with respect to the radial coordinate r. The equations 
of motion for an anisotropic fluid are [2JJ] 



47TP = L 



3 ( a + b ) Hr> 



4irp r 
4np t 



■ r 2 2 4 



fr 
2 

cot0 
2r 2 



T 



+ e~ 



, a 1 

2 + 4 + ^ 



b' 



I 

4 : 



T'fTT = 



(6) 
(7) 
(8) 
0) 



where p r and pt are the radial and tangential pressures respectively, p is density profile. This last 
quantity is very importnant in our astrophysical predictions. . Here if we use from the '"good"' 
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tetrads 22|, the out coming system becomes very complicated as the following: 



4vrp = L _ f -I^(2 -2e b + r 2 e b T - 2r b') - tlHl _ (l + e W sin^OLO) 



/ | /re- 
4 4r 2 



4vrp = -L + (2 - 2 e b + r 2 e 6 T - 2r a'flU) 



/ TT T' cos 7 =(CL2) 
/ T 6=(Q3) 

6[e fe r 2 /T + 2 /tt(1 + e b/2 sin 7 )(2 -2e b + r 2 e b T + 2r a')} - 4r f TT d'(l + e b/2 sin 7 ) 2 = (Q4) 
frr \ ~ 4 e a r T' - 6 2 (2 - 2 e b + r 2 e b T) - 2r a'(e a r T' + 6 2 ) + 4r6d'(l + e fe/2 sin7Xjl5) 
4 e a - 4e a e fe - e a r 2 a' 2 + 2 e a r b' + e a r a' (2 + r b') - 2r 2 e a a" - e b r 2 ab + eV 2 6 2 + 2e b r 2 'b\ = 



Here Jtt is the second derivatives of f(T) and overdots denote differentiation with respect to t, 
the Euler angles is j(r,t). We do not give any further discussion on this system, we will follow 
the simple, tetrads basis, given by ([!]). Note that even in this new tetrads basis, the teleparallel 
case is preserved by equation (112p . Thus even by this strange basis of tetrads, we can include the 
teleparallel case with linear T behavior. There is a vast family of exact solutions for this system, 
which has been investigated previously [20]. We focus only on the following possible solutions, 
which arise from both equations (|9ll2p : 

T' = =>• T = T , (16) 
fTT = => f(T) = C 2 + CiT , (17) 
T' = 0, f TT = => T = T , f(T) = /(T ); , (18) 

which always relapses into the particular case of teleparallel Theory, with f(T) a constant or a 
linear function. We adopt this linear teleparallel choice for our physical discussions about the 
possible explanation of the DM in the context of the torsion based gravity, f(T). In the next 
section, we will solve the above equations (J6]), ([7]) and ([8]) for the metric function a(r). In the 
language of the 3 + 1 decomposition of the metrics, determining a(r) is equivalent to finding the 

a(rj 

lapse (or redshift) function N(r) = e~^~ . 

IV. DARK MATTER PROBLEM IN f(T) GRAVITY 

The quasi global solution for ([3]) with assumption a(r) = —b{r) and by imposing the isotropicity 
in the pressure components p r = pt is the Schwarzschild-(A)dS presented by [3] 

e o(r) = e -b(r) = 1 _ 2° + ^i r 2_ 
r 3 
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Obviously such trial metric can not be successful to generating the rotation curve of the spiral 
galaxies. Indeed this classical solution leads the zero torsion T = 0. From physical intuition we 
know the DM problem must be comes from a non zero torsion, and specially from a variable one, 
T = T(r). Now we introduce an ansatz for solution and choose 

b(r) = c, (19) 

where c is an arbitrary constant 1 . Another choice is the polynomial form for the lapse function e a 



251 ]. but which such choices , independent from the origion which they come, the rotation curve 
fixed by a desirable linear form, and it seems that such choices are ad hoc and not physically 
acceptable. Further we choose another ansatz c\ = 1 and C2 = 0. It means f(T) = T. The 
model and the field equations still remain scale invariant. The main reason for choice of the 
metric function b as a constant goes back to the scale invariance of the system and further, we are 
interested to a lapse function N = e a which can explain the flat rotation curve. We discuss more 
on why we choose such a restricted gauge. Let us consider the following static form of the metric 

g^dx^dx" = h AB dx A dx B - $ 2 dE| (20) 

instead of ([3|), which is a four dimensional dual of the following renormalizable effective action, 
defined on the two dimensional induced metric 1iab(x, t), A, B = 0, 1, 

S = — [ d 2 xV Z h\f($)T + e$- A & A - U($j\ , e = constant 0(1). (21) 

167T J L J 

This action is a generalization of the action which is proposed in the Einstein gravity for some 

dilaton fields [24]. We assume that the free functions f, U are analytic in <3? in the limit of large <!>. 



Indeed this action is power-counting renormalizable 26]. This power-counting renormalizability is 
valid for any polynomial form of the interaction coupling f,U. Comparing ([3|). (|20p we find there 
exists a gauge freedon for choice the field $ = b(r). One trivial gauge is the constant gauge given 
by (|19p . Thus our ansatz can be interpreted asa gauge free term in the effective action. Now we 
back to the analytic investigation of the solutions. Without loss of any generalization,we assume 
the isotropic ansatz for the matter distribution 

Pt = Pv (22) 

The expression for rotation curves of galaxies is 

'r(e a )\i/2 



2e a 



(23) 



1 Here c does not represent the speed of light 
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where prime denotes differentiation with respect to radial coordinate r. This formula is the same 
as the Einstein gravity. We must clarify this point here. The path of the free particle can be 
obtained using the usual minimization method of the action for a free particle 

5 J ds = 0, ds = y / g^ iu x^x u drj. 

Here r\ is the affine parameter. Using this equation we obtain the following geodesic equation, 

x" + T£ a x u x a = 0. 

Here T„ a is the Levi-Civita connections, defined by the symmetric part of the general connection 
from the metricity equation. The asymmetric part of the connection has no portion in this 
geometrical equation. Thus if we mean by the geodesic equation, the non auto parallel motion, the 
same expression can be used. But using the auto parallel formalism is another story and we will 
not enter in it. It is easy to show that this geodesic equation is equivalent by the Euler-Lagrange 
equations, derived from the following point like Lagrangian for the test particle 

2C = e a i 2 -e b r 2 -r 2 (9 2 + sin 2 9<j> 2 ). 

The purely radial equation for test particle reads 

r 2 + C/(r) = 0, U{r) = e- h ^(e- a ^e 2 -%-l), h 2 =p 2 + -^-. (24) 

Here pe,Pip are the conjugate monenta of the corresponding coordinates, e = te a is the energy 
(local). Here we are assuming that there exist ZAMO (zero angular momentum ) observer, located 
at the spatial infinity r — > oo. We note that the tangential (rotation) velocity reads 



Mr) h 2 ^ ra /( r ) 



r 2 e 2 2 ' 

which coincides with (|23p . Making use of isotropic ansatz (|22p in ([7]) and ([8|), we obtain 



e aW = -^^(dr^ - C 2 )\ (25) 
16(-13 + 4e c ) 2V ; y ' 

where C\ and C2 are new constants of integration. A simple direct calculation of torsion using 
([5]) and with metric function given by (|25p shows that in this case T(r) 7^ (See figure 1). By 
substituting ([25]) in (f23|) we can plot the rotation (tangential) velocity as it is plotted in figure 
2. We used from a large set of data come from the local tests of the f(T) models based on the 



cosmographic description [29|. Figure 1 resembles the rotation curves for large (spiral) galaxies 
[30! . The scale where the velocity profile tv, attains is v v w 10~ 3 , which roughly corresponds 
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to r rs 300A;m/s, in reasonable agreement with the data. It is important to mention that the 
velocity profile shown in figure 2, is constructed basically on a phenomenological toy model. It 
is gratifying that our f(T) theory predicts a good velocity profile that was argued to be a good 
phenomenological fit to the data. 

io 7 - 

10 6 - 

io 5 - 

\n(T) 10 4 ' 

io 3 - 

10 2 - 

10 1 

0.2 0.4 0.6 0.8 1 

FIG. 1: Non-zero scalar torsion T for metric function given by (|25|) . 




V. ON CENTRAL DENSITY PROFILE p(r) IN f(T) MODEL 

From observational data we know that there exists a core with roughly constant density (mass 
density) in the galaxy. Many models have been proposed for this mass profile density. In brief 
these models are used in the numerical simulations, specially for ACDM: 

pW = m^w (26) 

Here pi represents the density of Universe at the collapse time, po is central density of the halo, r s 
is a characteristic radius for the halo and r c is the radius of the core 3l|. Further, very close to the 



center, this density profile is characterized by r n where n ~ — 1, which means we have a density 
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FIG. 2: Toy model rotation curve for large spiral galaxy based on our f(T) model. The units of the vertical 
axis must be read by as factor xlO-^ 2 . 



cusp. The observations, often favor n ~ 0, i.e. a constant-density core. Now we want to compare 
our estimation on the form of the p(r) from model of f(T) by metric function (]25f) . Substituting 
(1251) in (El) we obtain 



p(r) 



1 

16vr L 



c 2 + 



2ci{\-e 



(28) 



This is the exact mass profile density of our model of DM, which is comparable with the two 
recently proposed models of the density (|26|27|) (See the figure 3 for a comparison of the estimated 
model of us given in ([28]) and two models (|26|27p from astrophysical data given from [311] ). 



VI. CONCLUSION 



In this letter we obtained the rotation curve of the galaxies in the f(T) gravity. We proposed 
that the galaxy metric remains spherically symmetric and static. Then by solving the general 
equations of the metric components we obtained the lapse function N = e a as a function of the 
radial coordinate r. Our qualifying discussions shown a very good agreement between the rotation 
curves in this model and other curves obtained from the data. The scale where the velocity profile 
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login (p(»0) 




t 1 1 1 1 1 1 1 r~ 

0.2 0.4 0.6 0.8 

r(lOkpc) 



FIG. 3: Two astrophysical halo density profiles in comparison to our toy model. Solid is our prediction 
given by ([28)) , dot is the density of NFW model ([26| and dash the estimation of pseudo-isothermal sphere 
approximation introduced by (|27|) . The behavior of our model is very close to the estimations of the NFW 
model [311 ] . 

attens is v v ~ 10~ 3 , which roughly corresponds to r ~ 300 km/s, in reasonable agreement with 
the data. Further, the exact mass profile density of our model of DM, is in good agreement with 
the two models of the density. It proves that dark matter problem can be resolved as the effect of 
the torsion of the space time easily. 

It is well-known from ECSK theory that torsion couples to the spin of matter, 
hence one can infer to measure the torsion produced by any massive spinning object 
(including massive spiral and elliptical galaxies). However there has been assumptions 
about testing these theories, such as "all torsion gravity theories predict observation- 
ally negligible torsion in the solar system, since torsion (if it exists) couples onl y t o 



the intrinsic spin of elementary particles, not to rotational angular momentum" 32]. 
Mao et. al. have shown that Gravity Probe B is an ideal experiment for constraining 
several torsion based theories [321 ]. Although their analysis is based on the torsion field 
around a uniformly rotating spherical mass such as Earth, the task of constraining 
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torsion around massive galaxies is still open for exploration. 
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